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1. Introduction 

There is a profound diversity between quantum mechanical and classical correlations. 
Schrodinger [1, 2] coined the term "entanglement" to describe the pecuhar connection 
that can exist between quantum systems, that was first perceived by Einstein, Podolsky 
and Rosen [3] and has no analogue in classical physics. Entanglement is a resource in 
quantum information science [4, 5, 6] and is at the origin of many unique quantum 
phenomena and applications, such as supcrdense coding [7], teleportation [8] and 
quantum cryptographic schemes [9, 10, 11, 12]. 

Much progress has been made in developing a quantitative theory of entanglement 
[5, 6]. The bipartite entanglement between simple systems can be unambiguously defined 
in terms of the von Neumann entropy or the entanglement of formation [13, 4, 14]. On 
the other hand, an exhaustive characterization of multipartite entanglement is more 
elusive [5, 6] and different definitions [15, 16, 17, 18, 19] often do not agree with each 
other, essentially because they tend to capture different aspects of the phenomenon. 
More to this, a complete evaluation of global entanglement measures bears serious 
computational difficulties, because states endowed with large entanglement typically 
involve exponentially many coefficients. 

We proposed in [20] that multipartite entanglement shares many characteristic 
traits of complex systems and can therefore be analyzed in terms of the probability 
density function of an entanglement measure (say purity) of a subsystem over all 
(balanced) bipartitions of the total system [21]. A state has a large multipartite 
entanglement if its average bipartite entanglement is large. In addition, if the 
entanglement distribution has a small standard deviation, bipartite entanglement is 
essentially independent of the bipartition and can be considered as being fairly "shared" 
among the elementary constituents (qubits) of the system. Clearly, average and standard 
deviation are but the first two moments of a distribution function. A full characterization 
of the multipartite entanglement of a quantum state must therefore take into account 
higher moments and/or the whole distribution function, in particular if the latter is not 
bell shaped or is endowed with unusual and/or irregular features. 

The idea that complicated phenomena cannot be summarized in a single (or a few) 
number(s), but rather require a large number of measures (or even a whole function) 
is not novel in the context of complex systems [22] and even in the study of quantum 
entanglement [23] . In this article we shall pursue this idea even further and shall study 
the bipartite and multipartite entanglement of a system of qubits by making full use of 
the tools and techniques of classical statistical mechanics: we shall explore the features 
of a partition function, expressed in terms of the average purity of a subset of the 
qubits: this will be viewed as a cost function, that plays the role of the Haniiltonian. 
Interestingly, this approach brings to light the presence of frustration in the system [24] , 
highlighting the complexity inherent in the phenomenon of multipartite entanglement. 

This paper is organized as follows. We introduce notation and define maximally 
bipartite and maximally multipartite entangled states in Sec. 2. Multipartite 
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entanglement is characterized in terms of the distribution function of bipartite 
entanglement in Sec. 3. The statistical mechanical approach and the partition function 
are introduced in Sec. 4. The high temperature expansion and its first three cumulants 
are computed in Sec. 5. Section 6 contains our conclusions and an outlook. 



2. From bipartite to multipartite entanglement 
2.1. Bipartite purity 

We consider an ensemble S = {1,2, . . . ,n} oi n qubits in the Hilbert space Hs = (C^)®'^ 
and focus on pure states 

= Zk\k) , ZkeC, ^\zkf = l, (1) 



where k — {ki)i^s, with /cj e Z2 = {0, 1}, and 

\k)^(S)\ki)i, \ki)ieC^ {h\kj)^5ij. (2) 

ieS 

For the sake of simplicity, in this paper we shall focus on pure states of qubits and 
shall not discuss additional phenomena such as bound entanglement [25, 26]. Consider 
a bipartition {A, A) of the system, where A <Z S is a. subset of ha qubits and A = S\A 
its complement, with + = n. We set ua < tia with no loss of generality. The 
total Hilbert space factorizes into 1-Ls — 'Ha ® 'Ha-, with T-La = (S'ieA ' dimensions 
Na — 2"^ and Na — 2"^, respectively {NaNa — N). As a measure of the bipartite 
entanglement between the two subsets, we consider the purity of subsystem A 

7rA = trAPA> PA = tr^|V')(V'|, (3) 
trx being the partial trace over X — A or A. We notice that tt^ = 'ka ^-^d 

1/Na < tta < 1. (4) 
State (1) can be written according to the bipartition {A, A) as 

IV^) = X] Zk\kA)A ® \kA)A, (5) 

where kA = {ki)i^A and \1)a = (S)ieA ^ '^a- By plugging Eq. (5) into Eq. (3) we 
obtain 

PA= ^ ZkZl5kA,lA\^A){lA\ (6) 



and 



= 5Z ZkZk'ZiZi>5kA,l'Jk'^,lA^kA,lA^k'^,l'A' (7) 



which is a quartic function of the coefficients of the expansion (1). If, for example, 
the system is partitioned into two blocks of contiguous qubits {C,C), namely C ~ 
{1,2, . . .,nA), then 

TTC = ^ ^ Z(],^m)Z{V ,m)Z{l' ,m')Z{l,m')-i (8) 
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where {l,m) = (Zi, . . . , Zn^, mi, . . . , m„^) e Z^. 

2.2. Minimal bipartite purity 

For a given bipartition it is very easy to saturate the lower bound 1/Na of (7). For 
example, 

which represents a maximally bipartite entangled state 

\i^) = N2'^" E IOa^IOa, (10) 

yields Pa = I/^a and tta = '^/Na- In fact, the general minimizer is a maximally 
bipartite entangled state whose Schmidt basis is not the computational basis, namely, 

= ^~a'" E K,iuLi^ (11) 

where f/^"^, = {Ia\U'^\1'a) with f/^ a local unitary operator in I-La that transforms the 
computational bases into the Schmidt one, that is 

\^)-N-"" E U^\1)a®U^\1)a- (12) 

For this state, 7r^(| ■?/')) = I/Na- The information contained in a maximally bipartite 
entangled state with ua = ua is not locally accessible by party A or A, because their 
partial density matrices are maximally mixed, but rather is totally shared by them. 

2.3. Average purity and MMES 

Entanglement, in very few words, embodies the impossibility of factorizing a state of 
the total quantum system in terms of the states of its constituents. Most measures of 
bipartite entanglement (for pure states) exploit the fact that when a (pure) quantum 
state is entangled, its constituents do not have (pure) states of their own. This is, 
for instance, what we did in the previous section. We wish to generalize the above 
distinctive trait to the case of multipartite entanglement, by requiring that this feature 
be valid for all bipartitions. 

Let l'^) e Hs and consider the average purity [27, 29] 

^'Mkm-nu-(l) E (13) 

V ^/ \A\=nA 

where E denotes the expectation value, \A\ is the cardinality of A and the sum is over 
balanced bipartitions ua = [^/2], where [x] denotes the integer part of x. Since we are 
focusing on balanced bipartitions, and any bipartition can be brought into any other 
bipartition by applying a permutation of the qubits, the sum over balanced bipartitions 
in (13) is equivalent to a sum over the permutations of the qubits. The quantity ttme 
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measures the average bipartite entanglement over all possible balanced bipartitions and 
inherits the bounds (4) 

< < 1. (14) 

The average purity introduced in Eq. (13) is related to the average linear entropy 
Sl — Jl'ti i^ ~ ttme) [27] and extends ideas put forward in [18, 28]. 

A maximally multipartite entangled state (MMES) [29] \(f) is a minimizer of ttme, 

4t(l^)) = 4"\ (15) 
with = min{7rLt(|V)) I IV') e ^5, (m = 1}- 

The meaning of this definition is clear: the density matrix of each subsystem A C S oi 
a MESS is as mixed as possible (given the constraint that the total system is in a pure 
state), so that the information contained in a MMES is as distributed as possible. 

2.4- Perfect MMES and the symptoms of frustration 

For small values of n one can tackle the minimization problem (15) both analytically 
and numerically. For n — 2,3,5,6 the average purity saturates its minimum in (14): 
this means that purity is minimal for all balanced bipartitions. In this case we shall 
say that the MMES is perfect. 

For n = 2 (perfect) MMES are Bell states up to local unitary transformations, 
while for n = 3 they are equivalent to the GHZ states [30]. For n = 4 one numerically 
obtains eJ^^ = minvr^^ = 1/3 > 1/4 = 1/Na [29, 31, 32, 33]. For n = 5 and 6 one 
can find several examples of perfect MESS, some of which can be expressed in terms of 
binary strings of coefficients [zk = H in Eq. (1)]. 

The case n — 7 is still open, our best estimate being E^^^ ~ 0.13387 > 1/8 = 1/Na- 
Most interestingly, perfect MMES do not exist for n > 8 [27]. These findings are 
summarized in Table 1. This brings to light an intriguing feature of multipartite 
entanglement: we observed in Sec. 2.2 that it is always possible to saturate the lower 
bound in (4) 

TT?^ = 1/Na (16) 

for a given bipartition {A, A). However, in order to saturate the lower bound 

e!P^ = 1/Na (17) 

in Eq. (14), it must happen that (16) be valid for any bipartition in the average (13). 
As we have seen, this requirement can be satisfied only for very few "special" values 
of n. For all other values of n this is impossible: different bipartitions "compete" with 
each other, and the minimum Ej^^ of tt^^ is strictly larger than 1/Na- We view this 
"competition" among different bipartitions as a phenomenon of frustration: it is already 
present for n as small as 4 [24] . (Interestingly, an analogous phenomenon exists also for 
"Gaussian MMES", see [34].) 

This frustration is the main reason for the difficulties one encounters in minimizing 
ttme in (13). Notice that the dimension of Tis is N = 2"' and the number of partitions 
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Table 1. Perfect MMES for different number n of qubits. 



n 


perfect MMES 


2,3 


exist 


4 


do not exist 


5,6 


exist 


7 


? 


> 8 


do not exist 



scales like 2^. We therefore need to define a viable strategy for the characterization of 
MMES, when n > 8. 



3. Probability distribution of bipartite entanglement 

We now introduce the distribution function of purity over all bipartitions, pi'KAj-i that 
will induce a probability-density-function characterization of multipartite entanglement. 
For rather regular (i.e. bell-shaped) distributions the first few moments already 
yield a good characterization: in particular, the average will measure the amount 
of entanglement of the state when the bipartitions are varied, while the variance 
will quantify how uniformly is bipartite entanglement distributed among balanced 
bipartitions. 

The calculation of the properties of tt^ is particularly simple for an important class 
of states. Consider the set 



C 



|(^i,^2,...,^iv)eC^l5^ 

I k 



ZkV = 1 



(18) 



corresponding to normalized vectors m.'Ks- This set is left invariant under the natural 
action of the unitary group U{'Hs)- A typical state is obtained by sampling with respect 
to the action of U{Hs) on this set. Typical states have been extensively studied in the 
hterature [35, 36, 37, 38, 39, 40] and can be (efficiently) generated by a chaotic dynamics 
[41, 42]. 

For large A^, the tt^'s have a bell-shaped distribution over the bipartitions with 
mean and variance [21] 



N + l 



{{tta - l^f)o 



2{Nl 



1) 



(19) 
(20) 



(A^ + l)2(iV + 2)(iV + 3)' 

respectively, where the brackets (• • •)o denote the average with respect to the unitarily 
invariant measure over pure states 



TT 



N 



\zk\^ 1 JJd^^fedzfe, 

k ) k 



(21) 
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induced by the Haar measure over U{l-is) through the mapping — ^ Zj\j) — C/|'0o), 
for a given reference state \iPq) [38]. Here dz^d^jt = dx^d^fe, with — Ylez^ and 
yk — Im^fc, denotes the Lesbegue measure on C. 

Given a state |^) G "H^, the potential of multipartite entanglement has the following 
expression in terms of its Fourier coefficients Zi 

ttme^ A{k,k';l,l') ZkZk' zizi< , (22) 

k,k',l,l'&Z^ 

with a coupling function 
A{k,k';i,i') = I M J2 I ('^feA,i:4^^:4.'A'^fcA,'A'^fc:4''A + '^^A.^:4^*=^.'A4;5,i^4^,Z^) , (23) 

\ ^/ \A\=nA 

with riA — [n/2] (balanced bipartitions) . The result follows by plugging the expression 
(7) of 71 A into Eq. (13), and by symmetrizing under the exchange k ^ k' (or, equivalently, 
A -H- A). The coupling function A has the following expression (see Appendix A for 
details) 

A(A;, k'; I, I') = g{{k ® I) V {k' © I'), {k I') V {k' I)), (24) 

where 

g{a,b) ^ 5aAb,o 9{\a\,\b\), (25) 

with \a\ = ^jg^Oj, \b\ = JZi&s^i^ a Q) b = {ai + bi mod 2)jg5 is the XOR operation, 
aV b — {ui + bi — aibi)i^s the OR operation, a Ab — {aibi)i^s the AND operation and 

^, , 1 ( n\ ( n — s — t\ ( n — s — t 
2 yiA J [\ riA- s J y UA-t 

Using the definitions we notice the following symmetries of the coupling function: 

A{k,k';l,l') = A{k',k;l,l') 

A{k,k';l,l') = A{l,l'-k,k') . (27) 
A{k,k';l,l') = A{k',k;l',l) 

4. Pcirtition function 

In order to study the minimization problem, we will reformulate it in terms of classical 
statistical mechanics: in particular, the minimum Eq of ttme will be recovered in the 
zero temperature limit of a suitable classical system. 

The main quantity we are interested in is the average bipartite entanglement 
between balanced bipartitions, ttme in Eq. (13). This quantity will play the role of 
energy in the statistical mechanical approach. We therefore start by viewing ttme in Eq. 
(13) as a cost function (potential of multipartite entanglement) and write 

H{Z) = TTMEd^)), (28) 

where z are the Fourier coefficients of the expansion (1). We consider an ensemble 
{rrij} of M vectors (states), where rrij is the number of vectors with purity H = Cj. In 



(26) 
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the standard ensemble approach to statistical mechanics one seeks the distribution that 
maximizes the number of states Q = M\ / Y[j under the constraints that mj — M 
and J2j '^3^3 = ME. For M — )■ oo, the above optimization problem yields the canonical 
ensemble and its partition function 

= j <\^xc{z)e-^"^'^ ^ CN j d//H(C/)exp(-^E[tr^(tr^C/|V'o)(V'o|C/^)']), (29) 

where the expectation value E was introduced in Eq. (13) and the Lagrange multiplier 
/3, that plays the role of an inverse temperature, fixes the average value of purity 
E. In the first integral we have used the measure (21) and taken into account the 
normalization condition (18). In the last (base- independent) expression //h denotes 
the Haar measure over U{l-L), jV'o) is any given vector and the (unimportant) constant 
ctv is proportional to the ratio ;U(7(C^)//iH(W('H)) between the area of the {N — 1)- 
dimensional sphere (18) and the volume of the unitary group. In conclusion, the 
potential of multipartite entanglement can be now considered as the Hamiltonian of 
a classical statistical mechanical system. 



Comments 

In order to clarify the rationale behind our analysis, a few comments are necessary. 

i) Although our interest is focused on the microcanonical features of the system, 
namely on "isoentangled" manifolds [43], we find it convenient to define a canonical 
ensemble and a temperature. This makes the analysis easier to handle and is at the 
very foundations of statistical mechanics, when one discusses the equivalence in the 
description of large systems between the microcanonical ensemble (in which energy is 
fixed) and the canonical ensemble (in which temperature is fixed). 

ii) One can view the multipartite system as an ensemble for the collection of all 
balanced bipartitions. However, what makes the problem intricate and interesting is the 
fact that there is a nontrivial interaction among different bipartitions, which in general 
provokes frustration. 

iii) From a physical point of view, the measure of typical states is a uniform measure 
over the whole projective space. This would be consistent with ergodicity. However, our 
analysis is purely static and we are not considering the time evolution generated by the 
(purity) Hamiltonian. The relaxation to equilibrium, as well as its ergodic properties, 
deserve a deeper study and would probably uncover additional features with respect to 
the equilibrium situation. This aspect will be investigated in the future. 

iv) Temperature is a Lagrange multiplier for the optimization parameter. It is the 
variable that is naturally conjugate to i/, in exactly the same way as inverse temperature 
is conjugate to energy: /3 fixes, with an uncertainty that becomes smaller for a larger 
system, the level of the purity of the subset of vectors under consideration, and thus an 
isoentangled submanifold. The use of a temperature is a common expedient in problems 
that can be recast in terms of classical statistical mechanics. One can find examples of 
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this kind in the stochastic approach to optimization processes (for instance simulated 
anneahng) [44, 45]. 

4-2. Some limits 

Wc start by looking at some interesting limits and give a few preliminary remarks. 
For /3 — )■ 0, Eq. (29) clearly yields the distribution of the typical states (21). For 
P — )■ +00 (T — )■ O"*"), only those configurations that minimize the Hamiltonian survive, 
namely the MMES. There is a physically appealing interpretation even for negative 
temperatures: for /3 — oo (T — > 0^), those configurations are selected that maximize 
the Hamiltonian, that is separable (factorized and non-entangled) states. 

The energy distribution function at arbitrary (3 can be obtained from the partition 
function 

Z{(3) = I d/ic(^)e^''^(^) = /' dEe-^^ / d^c{z) 5{H{z) - E), (30) 

J J Eo J 

where E G [Eq, 1], Eq being the minimum of the spectrum of H and 5 the Dirac function. 
Incidentally, notice that Eqs. (14) and (19) yield 

hm < lim = 0, > 2-["/2l, /x^") = = — . 

n—^oo n— >-oo 2" -|- 1 

(31) 

The energy distribution function reads 

^^^^^ ^mJ "^^""^'^ ^^^^'^ ~ ^^^^ 

which, for /3 = 0, simply reads 

^ ^ / '^'"''^'^ ^^^^'^ ~ ^^^^ 
In Fig. 1 we show the probabihty density function Po{E) for n — A. As emphasized 
in Sec. 2.4, this is one of those cases in which frustration appears, as for n = 4 qubits 
one (numerically) finds E^"^^ = minTT^J^ = 1/3 > 1/4 = 1/Na [29, 31, 32, 33]. We 
clearly observe the asymmetry of the curve, denoting a positive value of the skewness. 
This deformation becomes less evident for larger values of n. As we will see, in the 
thermodynamic limit, n oo, Po{E) will become more and more symmetric and will 
tend to a Gaussian. 

Using Eqs. (30)-(33) we obtain the expression of the energy distribution function 
at arbitrary 13 in terms of its infinite temperature limit: 

= . (34) 

Sl^iEe-l>^P,(E) 

Notice that this equation is valid at fixed n. 

By multiplying and dividing the last equation by |/3|e^ and /3e^^°, respectively, and 
remembering that 



1 




Figure 1. (Color online) Probability density function Pq{E) for n ~ A. The 
distribution has been obtained from 5 x 10''^ typical states. The binning is 3 x 10^^ 
and the integral is equal to 1. 



we find 

P_oo(^) = 5(^-1), Poo{E) = 6{E-Eo). (36) 

These limits are the counterparts of those discussed for the partition function and are 
reflected in the asymptotic behaviour of the average energy as function of /3 

{H)p = ^/ ^l^c{z)He-P^ = £ dEEP^iE) = -^ln^(/3). (37) 

Indeed, 

{H)i3^-oo = 1, {H)i3^+oo = Eq. (38) 

More generally, the m-th cumulant of H reads 

am am— I 

-rt^] = i-r^^r.zi(i) = {-r-'^,{H),. (39) 

We find 



§^{H), = -nf[H] = -{H'), + {H)l ^ -4 ^ 0, (40) 

which is non-positive. In particular 

a' = al = k^^\h]. (41) 



The average energy is a non-increasing function of /3 and has at least one inflexion point 
as function of /3. Moreover 

- - (42) 
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From the qualitative behaviour of one can obtain information about the width of 
the distribution. For (3 — >■ +00 the curvature of {H)^ is positive and therefore is a 
decreasing function. 

From a quahtative point of view, one expects the behavior sketched in Fig. 2: 
for /3 — )■ 0"*" (T ^ +00), the distribution is beU-shaped (typical states); when 
P — )■ +00 (T ^ 0+) the distribution tends to become more concentrated around Eq. The 
energy distribution (34) at sufficiently high temperatures [how high will be discussed in 
Sec. 5.4, see Eq. (134)] can be obtained by observing that from Eq. (40) 



Pl3{E)r^Po{E + /3a^). (43) 

For larger values of (3 the left tail of the distribution starts "feeling" the wall at Eq. 
The value of Po(-^o) influences the behaviour of Pjj^E). In general, Pq{Eq) can vanish 
or not, yielding the behavior sketched in Fig. 2(a) and (b), respectively. One finds 

Pp{E) ~ ^{E - Eoye-^^^-^°\ (44) 

where r is the order of the first nonvanishing derivative of Po{E) at Eq. [Figs. 2(a), (b) 
display the case r = 0, 1, respectively] Notice that the only relic of Po{E) in (44) is r 
and Pj3^oo{E) yields the second equation in (36). Actually if r = 0, Eq. (44) yields a 
pure exponential converging to 

P+^{E)^5{E-Eo). (45) 

If r > 1 the probability for finite (3 has an initial polynomial increase but still converges 
to a Dirac 6 in £'0, corresponding to MMESs. The analysis for /3 — >■ —00 is analogous (we 
expand P/^iE) around E = 1, which is the maximum of H); it yields the first equation 
in (36). In this limit we obtain the separable states. 



5. High temperature expansion 

This section is devoted to the study of the cumulants of Po{E). This will enable us to 
look at some properties of the high temperature expansion of the distribution function 
of the potential of multipartite entanglement. We remind that for ^ — > one gets the 
typical states. 

The high temperature expansion originates from the Taylor series 

13=0 

-E^-&'[^]- (46) 

j=o 

The average energy reads 

(_fl\m-l 
m=l ^ ' 
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Figure 2. (Color online) Qualitative sketch of Eq. (34), at fixed Na, in arbitrary 
units. The energy density function is distributed around fi with standard deviation 
CT at /3 = (first bell-shaped curve at the right in both panels) and moves toward 
Eq when /? increases. In each panel, from right to left, (3 changes in constant steps. 
The probability density rigidly shifts with /3, for /3 < iV7/2-i°g2 3. Eq. (133) and 
following discussion. Note that both Eo,fi are 0{N^^^^). In (a) Po{Eq) ^ 0; in (b) 
PoiEo) = 0. 



~ {H)o -f3{iH- {H)of)^ + y ((^ - (^)o)'>o ' (47) 
while the free energy takes the form 
F(/3) = ilnZ(/3) 

~ ^ - {H)o + f ((^ - {H)o)\ - y ((^ - (^)o)^>o • (48) 

In the following three subsections we will evaluate the first three cumulants of the 
distribution for /3 = in order to characterize the high temperature expansion of the 
energy distribution function. 

5.1. First cumulant 

The joint probability density oi z = (zk) G associated to the measure of typical 
states (21) is 

Pm{zi,Z2,...,zn) = V IzfclM . (49) 




By integrating out N — M variables, one gets 

{N-iy. 



PNizi, Z2, . . . , Zm) 



{N-M -ly.TT 



M 




l<k<M 

for 1 < M < A^. In particular the probability density of an arbitrary element of z is 

P^(^i) = — (i-kilT^^. (51) 
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Since (e*^'^^^^ )o = 0, the only nonvanishing averages of the type 

\kex leY I 

with X,Y C Z2 , are obtained when the variables {zk} and {zi} are equal pair by 
pair, that is when the sets of indices are equal, X — Y . The nonvanishing correlation 
functions are given by 

Ik \ k 

n 1^^/"^' ) = / n i^.f '"^PAr(^i, . . . , ^iv) n dz,dz, 
\j=i / j 

{N-iy.uU^,^- ,,,, 

(53) 



A simple proof goes as follows. Extend the product to all N variables by letting 
some rrij vanish and consider the quantity, with ctj > 0, 

rr|^.|2m,e-E.a,|.,|A _ f rr(da;,-x7^)e-S^"^^^(7V-l)! f — e--(^-S*-^) 

= (7V-1)! / ^e--n>^m.K-i^)' (54) 



where 



J^{z)^ [ x"'e-^''dx, {Rez>0). (55) 

J'S.+ 

Now, we have Jo{z) = 1/z and Jrnl^^) = (— l)'"d™ Jo/d^; = m!/2;"*+^, and thus 

TT|^ |2rn,e-E.a.k.lA ^(TV-i)! / — e-'-TT . ^"^^'l' ^1 (56) 



J I 

By setting aj — a for all j we get 



2mj 

\ '0 



l^il 



= (^-l)!nK)!/^e-^— (67) 

which when mj = for all j reads 

e-«.(JV-l,./^^e-^^. (58) 

Therefore, 

/-i-r o \ i)!n,K)! 

W /o (^ + E,^.-i) 

and (53) follows. 
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The average energy at (3 — can be easily evaluated and is equal to the 
purity II defined in (19): 

Let us check the above result by direct computation, through (22). We get 

{H)o= A{ki,k2;li,l2){zkiZk^zi^zi^)Q. 



Now, 



and thus 



+ {{\zi\'^)o - '^{\zi\^\z2\^)o) hiM^kuhhiM 

ki,k2&2 

+ {{\z,W-2{\z,\^\z2W)Y,HKk;k,k), 



where the symmetry (27) was used. 

By using (24) and by setting A; = /ci © A;2, wc get 

^{h, k2] k,, k2) = ^(0' ^1 ® ^2) = J] g{K ^)^nY 9{K 0). 

ki,k2 k\,k2 k,k2 k 

Since Sk^o = 1 and Yjo<s<n^\k\,s = 1> by using (25) we can write 

5]^(A;,0) = J]^(|A;|,0) = J]^(|A;|,0) J] 5\k\,s^ Y 9{s,0)Y^\k\,s- 

k k k 0<s<n 0<s<n k 

The number of strings containing s ones is 



Y ^ii'i'' 



n 



and from (26) 



n 



9{s,0) 



1 I n\ I n 
2\sj [uA 

+ 



riA 
s 



whence 



Ua 

s 



-1 



s 



n, 



n — s 
riA- s 



n — s 
riA 



where Na = 2"-^ and A^^ = 2"^. Therefore, one gets 
Y Hki,k2;k,,k2) ^ NY9{k,0) 

k\,k2 k 



'-{Na + Na), 



N{Na + Na) 
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On the other hand, 

because ^(0, 0) = 1. Summing up, we get 

{H)o = 7V(7Va + Na){\zi\''\z2W + N {{\ziW - 2(|^i|2|^2r)o) , (71) 
and since [see Eq. (53)] 

(i..n.2p)o=i{i.:r)o=^^^^, (72) 

we obtain 

(^)o - (73) 

which equals the value (19) of the average purity /x*^"-*. In the thermodynamic limit, 
N ^ oo, with Na = Na = \fN 

(H). ~ ^. (74) 

5.2. Second cumulant 

The second cumulant is defined as 

a^ = Kf[H] = {H')^-{H)l (75) 

In order to evaluate this quantity we will use a diagrammatic technique based on the 
definition of the coupling function A and its properties [Eq. (27)]. We start considering 

(-f^^)o = E A(^i'^2;/i,/2)A(A;3, A;4;/3,/4)(^fc^2;fe2^fc3^fc4%^/2^/3%)o- (76) 

We must have {ki} = {Ij} as sets, that is 

k = kp(^i), P&Vi (77) 
with 1 < i < 4, where Vi is the permutation group of {1, 2, 3, 4}. Therefore, 
(if^)o = E E A(/ci,A;2; A;p(i),A;p(2))A(/c3, A;4;/cp(3),/cp(4))((|2;fei|^|2;fe2n2;jk3H2;fc^ 

where 

|2m|„ |2n|„ |2s|„ |2t\\ /|„ |2m|„ |2n|„ |2s|„ |2t\ ('7n^ 

\Z2\ \Zz\ \Zi\ ))o = — , ■ P2 ^4 )o- (79) 

The above normalization takes into account the fact that if ki — kj for some i j 
the sum over the permutation group V4 overcounts the number of different terms. For 
example, if ki — k2 and different from the others, we get m — 2, n — 0, s — t — 1, and 
there is a factor 1/2!, while, if ki — k2 — k^ ^ k^, we get m — 3, n — s — 0, t — 1, and 
there is a factor 1/3!. 

Since m + n + s + t — 4, from Eq. (53) we observe that 

((l^il \Z2\ \zs\ \z,\ ))o - - + i)(Ar + 2){N + 3) 
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is independent oi k E Z|". Therefore, 

{H')o = (i^iri^2ri^3ri^4r)o e [^(i) h2),p(3) p(4)], (si) 

with the notation 

[p(l) p(2),p(3) p(4)] = ^ A(A;i, k2; A;p(2))A(A;3, A;4; /c^o), ^^^(4)). (82) 

Note that by the symmetries (27) of A, we can swap p(l) <H- p(2) or p(3) •H- p(4), as 
well as 1 <H- 2 or 3 ^-j- 4, so that 

[w x,y z] = [x w,y z] = [w x,z y] = [x w, z y], 

[w x,y z\ = [y x,w z], if w,|/e{l,2}, or w,7/e{3,4}. (83) 

Using these symmetries we can give a graphical representation of the quantity in Eq. 
(82). Let us consider Fig. 3a. Each vertex represent a pair {ki,kj) in the summation. 
The edges between vertices and the loops on the same vertex fix the value oip{i). For 
instance, a double loop on {ki, ^2) and {k^, k^) (see Fig. 3b) yields 

[2 1,3 4] = J]A(fci, k2; k2, ki)A{k3, k^; k^, k^). (84) 

k 

Each vertex has order 4 with two incoming and two outgoing edges. Each graph 
is oriented. However, for simplicity, in the graphs of Fig. 3 wc have not indicated 
the orientations, since in this case, as it is easy to see, they do not yield different 
contributions. As wc shall sec, this will not be the case for higher cumulants, where 
graphs with different orientations represent noncquivalcnt contributions. 

We start considering graphs with no links between the left and right pairs, see 
Fig. 3b. The sum of this class of graphs is 

[0 - link] = [1 2, 3 4] + [1 2, 4 3] + [2 1, 3 4] + [2 1, 4 3] = 4 [1 2, 3 4]. (85) 

We have 

[1 2, 3 4] = 5^A(A;i, k^] ki, k2)A{k3, k^; k^, k^) = ^ ^(0, k^ ® k2) ^(0, k^ © k^). (86) 

k k 

By setting l2 — ki® k2 and I3 — ki ® k^, we get 

[12,34] = ^^ <7(0, /2) <7(0, h) =N'Y, 9(0, 12)9(0, k) = N^{Y1 ^(0' ^)) (8^) 

ki,k4 hjh hjh k 

Therefore, by using (68), we get 

[0 - link] = 4 [1 2, 3 4] = N'^{Na + N^f. (88) 

Let us now consider the graphs with two links between left and right pairs in Fig. 3c. 
The sum of this class of graphs is 

[2 - link] = [1 3, 2 4] + [1 3,4 2] + [3 1,2 4] + [3 1, 4 2] 

+ [2 3, 1 4] + [2 3,4 1] + [3 2, 1 4] + [3 2, 4 1] 

+ [1 4, 2 3] + [1 4, 3 2] + [4 1,2 3] + [4 1, 3 2] 

+ [2 4, 1 3] + [2 4,3 1] + [4 2,13] + [4 2, 3 1] 

= 16 [13, 2 4]. (89) 
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(ki,k2) 


(k3,k4) 


• 


• 




(a) 



(k],k2) 


{h,k4) 




(c) 




{ki,k2) 


^\ {ki,k4) 




(d) 



Figure 3. (Color online) Graphs used for the evaluation of the second cumulant (75). 
(a): empty diagram, (b): graph with no links between the left and right pairs, (c): 
graph with two links between the left and right pairs, (d): graphs with four links 
between the left and right pairs. 



One gets 

[1 3, 2 4] A(A;i, k-,- h, h)Aih, h; k^, k^) 

k 

= '^g(k2® k3, {ki © k^) V {ki © /cs)) g {k2 © /cs, (/c2 © k^) V (/cg © k^) 

k 

(90) 

By setting /i = fci © /cs, ^2 = ^2 © ^3, and 1^ = k^® k^ we get 

[13,2 4] g [k, k V (/i © /2)) g (/2, {h © U) V /4) 

= N J2 aikjiyk) gihhyU), (9i) 

Zl ,l2 ,1a 

where we have used the (easy to prove) useful relation 

k V ill © h) = hy k- (92) 

We get 

k A (/i V k) = {k A k) V k, (93) 

so that the constraint of the function g, k A {k V ^2) = 0, implies that k = 0. Therefore, 
by using (68), we obtain 

[2 - link] = 16 [1 3,2 4] = 16A^ ^ g{k, 0)) ^ = m{NA + N^f. (94) 
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The contribution of the graphs with four hnks between left and right pairs (see 
Fig. 3d) has the form 

[4 - hnk] = [3 4, 1 2] + [4 3, 2 1] + [4 3, 1 2] + [3 4, 2 1] = 4[3 4, 1 2]. (95) 

We have 

[3 4, 1 2] = 5^A(A;i, /ca; h, h)A{k3, h; h, = 5^A(A;i, /cs; h, k^f 

k k 

= ^c/(|(fci © ks) V (fcs © k^), {ki © ki) V (/cs © /cs)) • (96) 

k 

By setting l-^^ — ki® ks, U — ki® k^, and Z2 = ® ^2 ® ^3 ® ^4, we get 

[3 4,12] = J] J] 9{hy hf - NY,9{hM)\ (97) 

where we used the relation (92) and the constraint I2 — implied by (ZiVZ2)A(Z4VZ2) = 0. 
Therefore, we get 

[4-link] =4[3 4,12] =iV/2(iV), (98) 

where 

/2(iV)=4 ^ g{k,lf. (99) 

Notice that if 

d^\AnB\^\BnA\e[o,nA]. (100) 

is the distance between bipartitions {A, A) and {B,B), defined as the number of qubits 
belonging to A and not to B, then 



-1 



f,{N) = ^y^ j E yd) [dj t^''^'"' + 4-^'^/^-'^)] . (101) 

See Appendix B. Summing up, we obtain 

[p{l) p(2),p(3) p(4)] = [0 - link] + [2 - link] + [4 - link] 

= 4[1 2,3 4] + 16[1 3,2 4]+ 4[3 4, 1 2] 

^ N{N + 4){NA + NAf + Nf2{N). (102) 

Therefore, 

^^^0- (Ar+l)(7V + 2)(7V + 3) ' ^^'^^ 



and 



.2 (iV+l)/2(iV)-2(iV^ + iV;i; 



(iV+l)2(iV + 2)(iV + 3) ■ 
We have checked that the above analytic expression of the second cumulant, with /2 
given by Eq. (99), agrees very well (within 1% up to n = 8) with the numerical estimates 
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Figure 4. (Color online) Difference between the analytic value of second cumulant, 
computed according to Eq. (104) with /2 given by (99), and its asymptotic limit, 
obtained by substituting (105) into Eq. (104). The approximation is valid within a few 
percent even for n — 3. 



based on the probability density function (obtained by sampling 5 x lO'' typical states 
for each value of n). 

Finally, one proves that (see Appendix B), in the limit — ?■ oo, 

/2(A^) ~ 3y2iV", (105) 

with 

a = logaS - 1 ~ 0.5850. (106) 
Therefore, for — t- oo we have 

^^-4^-]^-0(iV— ). (107) 
Incidentally, note that 

-2 



a 




-2 
'^ind 




J2 [(^^^b)o - Mo Mo] , (108) 

A,B 

SO that, if the bipartitions were independent, we would have obtained 

E [«)« - <-4] = (,,; J « 

Thus, the result in Eq. (107) detects an interference among different bipartitions. We 
stress that the asymptotic estimate is very accurate even for small values of n. In Fig. 
4 we plot the difference between the analytic value of second cumulant, obtained using 
Eq. (104) with /2 given by Eq. (99), and its asymptotic limit, obtained by substituting 
(105) into Eq. (104). We notice an oscillatory behavior: the asymptotic expression 
systematically overestimates (underestimates) the second cumulant for even (odd) values 
of n. On the other hand, the approximation is very good even for small values of n. 
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5.3. Third cumulant 

The third cumulant is defined as 

k^^\h] = {{H - {H)of), = {H\ - 3{H\{H), + 2{H)l. (110) 
In analogy with the evaluation of the second cumulant we have 

{H% = (i^irk2ri^3ri^4rk5ri^6r)o e b(i) p(2),p(3) p(4),p(5) p(6)], 



fiir 



with 



b(l) p(2),p(3) p(4),p(5) p(6)] = ^ A{h,k2;kp^i),kp(^2)) 



X A{k3,k4; /cp(3), kp^4))A{k5, ke] /cp(5), A;p(6)) 

(112) 

and 7^6 the permutation group of {1, 2, 3, 4, 5, 6}. Prom Eq. (53) we easily obtain 



I |2| |2| |2| |2| |2| |2\ _ 
|^l| p2| |^3| F4| |^5| pel /O — 



N{N + 1){N + 2){N + 3){N + 4)(iV + 5) ' 

(113) 

We start by considering connected graphs with three ears. A representative of this 
equivalence class is depicted in Fig. 5a. We have 

[1 6, 3 2, 5 4] = J2 ^(^1' ^2; ^1, h)Hk3, kr, k3, k2)A{k5, h, k^) 

k 

= ^g(k2® ke, {ki e ke) V {ki /ca)) g(k2 ki, {k2 k^) V {k^ ® ki)^ 

k 

X g(k4® ke, {k^ ® /cg) V (/cg © /ce)) 

= ^ g{0, ki ® A;2) ^(0, k2 ® /cs) ^(0, /ca ® /cs) 

= iv(5:,(0..,)^iv(^^!^. (114, 

k 

where the constraint in the definition of the function g has implied k2 — k4 — ke and we 
used Eq. (68). The degeneracy of this class of graphs is 128. 

We now consider the class of connected graphs with two ears represented in Fig. 5b. 
We obtain 

[1 3, 2 5, 4 6] = ^ A{ki, k2; ki, ks)A{k3, k^, /cs, k5)A{k^, ke, k^, ke) 

k 

= 9{ki®k2,Q)g{k2®k4,0)g{k4®k6,0) 

ki,k2,k4,ke, 



N 

k 
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Figure 6. (Color online) Oriented connected graphs used for the evaluation of the 
third cumulant. (a): same internal and external orientations of the edges (nonvanishing 
"current"), (b): opposite internal and external orientations of the edges (vanishing 
"current"). 
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where we have imposed /c2 = and = k^. The degeneracy of the class is 192. 

The final class of non-oriented connected graphs is represented in Fig. 5c. Its 
explicit calculation yields 

[1 6, 2 5, 3 4] = ^ A(A;i, fcs; ki, fc6)A(/c3, k^; k2, kr^)A{k^, ke; fcg, k^) 



^g{k2® ke, {ki © fee) V {ki © k2 



X g[{k2 © fcs) V {ki © k^), {k^ © kr,) V (^2 © k^) 
X ^(^(A;3 © /cg) V {k4 © A:6), (A;4 © /^s) V (fcg © fee) 
= ^ ^(0, fci © k2)g({k2 © /cs) V {ki © feg), (fcg © k^) V (/ca © 



2 

4; ^ 



A:i,...,fc5 

= g{Q,ki®k2)g{k2®k^,k2®kif 

fcl,...,fc5 

= iV 5^ ,(0, k) Y: 9{hMf = N^^^±^ Y hf 

k lij.2 liM 

= N^^^^f2{N), (116) 

where we have used the constraint ^2 = ^6 and the function f2{N) defined in Eq. (99). 
The degeneracy of this graph is 192. 

In order to take into account the contribution of connected graphs with no ears, it 
is necessary to consider two different classes of oriented graphs whose representatives 
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are shown in Figs. 6a and 6b, respectively. For the first class (Fig. 6a, nonvanishing 
"current") we have 

[5 6, 1 2, 3 4] = ^ A(A;i, ki] h, ke)A{k3, k^; k^, /cs) A(A;5, k^; k^, k^) 

k 

= XI ^((^1 ® ^5) V {k2 e ke), (ki e ke) V {k2 k^fj 

k 

X g(^{ki © ks) V {k2 © k^), {k2 © /cs) V {ki © k^)^ 
X g(^{k3 © /cs) V (A;4 © ke), {k^ © /cg) V (/cg © /ce)) 

= gikl © fcs, ^2 © fcs) 5(^1 © ^3, k2 © /Cg) 

X g{k3 © /cg, /ci © /C2 © /ca ® ^5) 

Y giki,k2) giki®k3,k2®k3) g{k3,ki®k2®k3) 

J2 9{h,k2®h) g{k2,h®h) 9{h,ki®k2) 

k^MM 

^Nf^\N) (117) 
where we have used the constraints ki — k\®k2® kz and /ce = ^1 © ^2 © ^5 and defined 

fz\N)= 9{.kuk2®k3) g{k2,k^®k3) g{k3M®k2). (118) 

The degeneracy of this graph is 16. An analogous calculation can be carried out for the 
second class of oriented graphs (Fig. 6b, vanishing "current"). We obtain 

[3 6, 5 2, 1 4] = Y Hki. k2; ka, ke)A{k3, k4; k2)A(k5, ke; ki^k^) 

k 

Y 9{kiM)g{kMg{ki®k3M)^Nff\N) (119) 

with 

/f(^)= E 9{kiM)9{k3M)9{k,®k3M)- (120) 

In this case, the degeneracy is 64. 

The contribution of disconnected graphs (Fig. 7) can be computed by considering 
the results obtained for the first and second cumulant. For the class of graphs 
represented in Fig. 7a we have 

[1 2,3 4,5 6] ^7V3 (^^ + ^^)' , (121) 

8 

with degeneracy 8. In the case of the graph in Fig. 7b the result is 

11 3,2 4,5 61=iv(^^i±i!^iv(^5l+M.^.(^^l±^, (122) 
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with degeneracy 96. Finally from the disconnected graphs with two loops (Fig. 7c) we 
obtain 

[3 4, 1 2, 5 6] = N'i^^^^±Mf^^N), (123) 

o 

with degeneracy 24. In conclusion, we find 



{N + l)(iV + 2){N + 3)(iV + 4)(iV + 5) 
X (l6/f (^) + 64/f (iV) + 3{N + 8){Na + Na)MN) 

+ (Na + Na)\N^ + 12N + 40)) (124) 

and, therefore, 

^(3)ro-i = I 

° ^ J (AT + 1)3 (AT + 2) (AT + 3) (AT + 4) (AT + 5) 

X (l6(Ar + irfi'\N) + 6A{N + iyft\N) 

- 36{N + 1){Na + Ar^)/2(iV) - 8(Ar^ + Ar^)'(A^ - 5)) . (125) 

We have checked that the above analytic expression of the third cumulant, with /2, /g^^ 
and /g^^ given by (99), (120) and (118), respectively, agrees very well (less than 1% for 
n = 1-^7 and a few % for n = 8) with the numerical estimates based on the probability 
density function (obtained by sampling 5 x 10"^ typical states for each value of n). 
Finally, in the limit A^ — )■ oo, one can prove that (see Appendix C) 

fP{N)r.cN'-\ (126) 

with c ~ 1.05385 and 7 ~ 4.1583 given by (C.30) and (C.31), and that (see Appendix 
D) 

fi'\N) ~ AT", (127) 

with a ~ 0.5850 given by (106). Therefore, by recalling the asymptotic expression for 
/2(A^) (105), we find that the graph in Fig. 6b dominates over that in Fig. 6a and 

[H] ~ 64 c AT-T ~ 67.443 N'^-^^^^. (128) 

In Fig. 8 we plot the difference between the analytic value of third cumulant and 
its asymptotic limit, obtained by substituting (105), (126), (127) into Eq. (125). We 
again observe an oscillating behaviour. The approximation is good for n > 4. 



5.4- Gaussian approximation 

We can now summarize the results obtained for the first three cumulants and try to 
get a broader picture. Equations (73), (104) and (125) arc all exact. Their asymptotic 
expansions for large N, are given in Eqs. (74), (107) and (128). By plugging these results 
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Figure 8. (Color online) Difference between the analytic value of third cumulant, 
computed according to Eq. (125) with /a, ^ and f^^^ given by (99), (120), (120), 
and their asymptotic limits, obtained by substituting (105), (126), (127) into Eq. (125). 
The approximation is valid within a few percent for n > 4. 



into Eqs. (47) and (48) we obtain the asymptotic expressions of the average energy 
- 2 4^ ^67^ 



~ _/y0.5 ' _/Y2.415 ' 2 iV*^-1583 

and the free energy 

™)^l^_^ + ^^_^64c 
^' ' (3 ^/N 2 A^3-a 6 N'T 

InZ(O) 2 /3 4.243 /J^ 67.443 , , 

(5 2 Ar2.4i5 g ^^4.1583 V 

where 

c - 1.054, a = log2 3 - 1 ~ 0.5850, 7 = 4.1583. (131) 

See Eqs. (C.30) and (C.31). 

If is large enough and the first two cumulants at /3 = suffice, the energy 
distribution (33) can be taken to be Gaussian 

where /i and a are given in (74) and (107), respectively. The energy distribution at 
arbitrary temperature is then [see Eq. (43)] 

This is valid for relatively small /3: 

/X - /3cr' - (T > ^ /3<^/a2^iV^/2-iog,3_ ^^34) 
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Up to this value the probabihty density rigidly shifts with as is apparent in Fig. 2, 
which was obtained by numerically solving Eq. (34). 

5.5. A few comments 

The behavior of the cumulants derived in this section is very peculiar. The second 
and third cumulant follow a nontrivial power dependence, with trascendental exponents 
[see Eqs. (129)]). Interestingly, close scrutiny of the calculation in Sec. 5.3 shows also 
that 3 — a, the exponent that governs the A^-dependence of cr^, is found in a class of 
(nondominant) graphs that appear in the evaluation of the third cumulant: the exponent 
5 — a stems from the graph in Fig. 6a (the dominant exponent 7 stemming from the 
graph in Fig. 6b). This might suggest a possible recursion of the exponent a at all orders 
in the cumulant expansion. At this stage, we are unable to say if at higher orders the 
dominant graph for in Fig. 6b cancels, yielding a series in N'^^"'^ with q a function 
of a. 

It would be important to go beyond the Gaussian approximation in order to evaluate 
the behaviour of the left tail of the probability density function, close to ttme = — -Eq. 
See Figs. 1 and 2. This would give us some precious information about the features of 
MMES and the very structure of entanglement frustration [24] . In particular, it would 
be interesting to understand the role played by the interference among the bipartitions, 
in connection with the appearance of frustration in MMESs. See for instance the 
asymptotic behavior of the second cumulant in Eqs. (107)-(109) and the short discussion 
that follows. Additional investigation is necessary in order to elucidate these intriguing 
issues. 

6. Concluding renicirks and outlook 

We have built a statistical mechanical approach to multipartite entanglement, by 
introducing a partition function in order to tackle a complex optimization problem, 
whose solutions are the maximally multipartite entangled states, that appear as minimal 
energy configurations. 

The scheme adopted here is general. In classical statistical mechanics, temperature 
is used to fix the energy to a given value in the thermodynamic limit. Analogously, the 
fictitious temperature introduced here localizes the measure on a set of states whose 
entanglement (energy) is fixed, and can be larger or smaller than the entanglement 
associated to typical states. 

Remarkably, a strategy like the one adopted in this article, when applied to the 
simpler case of bipartite entanglement (at a fixed bipartition) [47] brings to light an 
involved landscape of phase transitions for the purity. Clearly, the multipartite version 
of the problem is much more involved, as the picture that emerges is complex and 
unearths a remarkable interplay between multipartite entanglement and frustration. 
It would therefore be of great interest to understand whether the phase transition that 
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occurs in the bipartite situation, when there is no average over the bipartitions, survives 
and has a counterpart in the multipartite scenario. This possibihty will be explored in 
the future. 

One important property that we have not investigated here and that is often used 
to characterize multipartite entanglement is the so-called monogamy of entanglement 
[15, 48], that essentially states that entanglement cannot be freely shared among 
the parties. Interestingly, although monogamy is a typical property of multipartite 
entanglement, it is expressed in terms of a bound on a sum of bipartite entanglement 
measures. This is reminiscent of the approach taken in this paper. The curious fact that 
bipartite sharing of entanglement is bounded might have interesting consequences in the 
present context. It would be worth understanding whether monogamy of entanglement 
generates frustration. 

Finally, we think that the characterization of multipartite entanglement proposed 
here can be important for the analysis of the entanglement features of many-body 
systems, such as spin systems and systems close to criticality. 
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Appendix A. 

We derive here the expression (24) of the couphng function. See [46] . We start from the 
definition (23), that can be rewritten as 

A{k, k'; I, I') = \ (^"^^ (A(fc, /, [n/2]) + K{k\ k; I, I'; [n/2])) ,(A.l) 

where 

A{k,k';l,l';nA) = ^ hj,,i'A^k'^,iJkA,iJk'^,i'^- (A.2) 

\A\=nA 

Let us fix a quadruple of binary strings {k, k', I, I') and a dimension ua- See figure Al. 
A bipartition {A, A), with \ A\ — ua yields a nonvanishing contribution to the sum (A.2) 
when 

that is when 

^A = I'a^ ^'a ~ f^A = Ia^ ^a = I'a^ (A-4) 

where we recall that kA = Ia means that the substrings of k and / arc equal, namely 
ki = li for all i E A. By noting that two bits ki and are equal when ki® h = 0, the 
above condition can be rephrased as 

kA®l'A = 0> k'j^®lA = 0, kA®lA = 0, k'A®l'A = 0, (A.5) 
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0000 


000000 


000 


00000 


0000 


000000 


111 


A, 


So 





Figure Al. (Color online) Graphic representation of the combination of binary strings 
contributing to sum (A. 2). 



that is 

{kA © I'a) V {k'^ © Ia) = 0, {kA © h) V {k'^ © I'a) = 0. (A.6) 

Summarizing, a bipartition {A, A) yields a nonvanishing contribution to (23) if and only 
if the following substrings are zero 

aA = and 6^ = 0, (A. 7) 

where 

a = {k®l)y {k' and 6 = (A; © /') V (A;' © /). (A.8) 

Note that equation (A. 7) implies that 

aA6 = 0, (A.9) 

since (a A 6)a = A = and (a A b)^ = A 6^ = 0. On the other hand, the substrings 
qa and Ba are totally free, whence 

\a\ = \aA\ < \A\ = ua, \b\ = {b^l < \A\ = n^. (A.IO) 

It is easy to see that (A.9) and (A.IO) are also sufficient conditions for the existence of 
a partition {A, A) that satisfies (A. 7). 

In conclusion, A(fc, k'; I, I'; Ha) 7^ when 

aA6 = 0, with |a| < n^, |&| < n^. (A.U) 

Therefore, 

A{k,k';l,l';nA) = 5aAb (|6|)#(fc,fc',/,r), (A.12) 

where xg is the characteristic function of set G and k', I, I') the number of terms 
in the sum in (A.l) that contribute to the function A. 
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Note that, since by (A. 11) the strings a and b cannot be both 1 at the same position, 
the set S is partitioned into three disjoint subsets (see figure Al) 

S ^ SoUAiUBi, (A.13) 

where 

^0 = {i e S\ai = 0, bi = 0}, 

Ai = {ie S\ai = l, bi = 0}, 

Bi = {ie S\ai = 0, bi = l}. (A.14) 

Obviously, 1^41 1 = \a\ < tia and = \b\ < n^- 

The number of terms ^{k, k' ,1,1') is given by the number of bipartitions {A, A) 
with \A\ = riA such that 

Aid A and Bi C A. (A. 15) 

Since ^4 fl ^4 = 0, parties A and A contend for Sq, namely 

A^ AiiJ{SQnA) and A^ B^iJ {Sq^ A). (A.16) 

Thus, the number of bipartition is the number of ways of picking |A\y4i| unordered 
outcomes from |S'o| possibilities. Since |^\^i| = |A| — |Ai| = tia — \a\ and |S'o| = 
IS"! — \A-i\ — \Bi\ —n—\a\ — \h\, one gets 

#(A;,^',M')= f (A.17) 
y ua- \a\ J 

Substituting Eq. (A.17) into Eq. (A. 12) and by defining the binomial coefficient to 
be identically zero when its arguments arc negative, wc notice that the characteristic 
functions in (A. 12) yield always one, and obtain Eq. (24). 

Appendix B. 

Wc derive here the asymptotic (for large A^) behavior of the function f2{N) defined in 
Eq. (105). Let us define the distance between bipartitions {A, A) and {B,B) as the 
number of qubits belonging to A and not to B 

d^\AnB\^\BnA\ e[0,nA]- (B.l) 
The number of pairs of bipartitions at a distance d is 

"-(:)(:)("/)■ 

Therefore the sum over the bipartitions can be rewritten as a sum over d 

\A\,\B\=nA \ ^ / d=0 \ / \ / 
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Let us consider for instance 



N 

-2 

N[nA 



-2 



fc \A\,\B\=nA 

+ 5(A;; i)5(A;; + S{k; A)S{k] B)^ 
J2 E l{s{k;A)6{k;B) + S{k;A)S{k;B)), (B.4) 



where 



S{k, k', I, I'; A) = dk^,i'jk'^^ijk^,ijk'^^i'_ (B.5) 
Let us start by showing that 

depends only on d — \Ar\ B\. When A — B, i.e. d = 0, 

A) = J] 5{k- Af = J2 S{k; A)= J2 ^'^A,vMu5kA,iJy,.v, = N\ (B.7) 

while, when d 7^ 0, we get 

h{A,B)=(^^J = 4-^N'. (B.8) 
Therefore, we get 

E S{k;A)5ik;B)^ J] hiA,B) 

k \A\,\B\=nA \A\,\B\=nA 

= e--*^-(:)e-^ (";)(";). (B.^) 

0<d<nA \ / 0<d<nA \ / \ / 

Analogously we find 

k \A\,\B\=nA \ J V,<d<nA \ J \ J 

Putting together Eqs. (B.9)-(B.10) we get 



k,l 

-1 



/ 5 



riA 



" > E ";)(";) (4"'" +4^) 

0<d<nA\ 
1 / 



= 2 ( „" j E (^7 I ( ) 2"" [4"""'' + 4-l»'''-'"] . (B.U) 
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We notice that the terms in the summation strongly depend on the ratio 2d/n. In 
the hmit n — )■ +00 only the terms with d = n/6 and d = n/3 give a significant 
contribution to the summation (see Fig. Bl). Let us consider the case of even n (in 
the thermodynamic limit the result for an odd number of qubits is the same) 

nA = nA = n/2. (B.12) 

By Stirling's approximation n! ~ (n/e)"'A/27m (for n large) and by defining the new 
variable 

2d 



X 



n 



after a straightforward calculation we obtain 

/ 7m „ 1 




2-"E 



exp {nS{x)}2 



n/2 



(B.13) 



2 nnx (1 — x) 

d 

1 

~ a/27™ / dx r 

Jo 27rx(l - X) 

X [exp{n[S(x) -X In 2]} + 2"" exp {n[5(a;) +x In 2]}] , (B.14) 

where 

S'(x) = — X Inx — (1 — x) ln(l — x) (B.15) 

is the Shannon entropy. Using the saddle point approximation in the integrand we get 

9 

4:71 



/2(A^) ~ V2^ I dx 
X exp < n S 



-\n2+-S" 
3 2 



1, 



X 



+ 
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v^Tm / dx— exp < n In - I 
Jo 47r ^1 2/ 



X 



cxp <^ ^"'^ I 



2^ f 9/ 2^'' 

+ exp <; -^^1 2; - - 



9 /3 



27r V2 



j j dxex.pl ~-xH ^3V2N'', 



where 



(B.16) 
(B.17) 



a = log2 3 - 1 ~ 0.584963. 
This is the asymptotic expression (105) used in Eq. (107). 

Appendix C. 

We evaluate here the asymptotic behavior of the function fl['\N) defined in (120). By 
using the definition (25), (120) can be written 

ft\N)= Yl 9ikiM)g{hM)9{h®hM) 

— <^feiAfc2,o5fc2Afc3,0^fe2A(feiefc3),0 

X ^(|A;i|, |A;2|) ^(|A;3|, |A;2|) g{\ki © \k2\) 

= XI /(«0,S1,S2,S3)^(S1,S2) ^(S3,S2) ^(Sl + S3 - 2so,S2), 



S0,Sl,S2,S3 



(C.l) 



where 



(C.2) 

kiMM 

and we have used 



/ci A A;2 = 0, A;2 A A;3 = =^ k2 ^ {h © h) = 0, 
\ki ® /csl = \ki\ + Iksl - 2\ki A A;3|. 

It is straightforward to count the number of terms in (C.2) and obtain 



(C.3) 
(C.4) 



f{so, 81,82,83) 



n \ / n - S2 \ I si \ ln-S2-Si 
S2J \ Si J \^oJ \ S3- So 



S2\so\{si - so)\{s3 - soy.{n - S2 - Si - S3 + So)! 
By substituting si — )■ si + sq and S3 — )■ S3 + sq in Eq. (C.l) we obtain 

n! 



. (C.5) 



S0,Sl,S2,S3 



so!si!s2!s3!(n - si - S2 - S3 - sq)! 
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X ^(si + So, S2) 9{s3 + So, S2) g{si + S3, S2) 
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E 



n 

So,Si,S2,S3 



where 



n 

11,12, ...,ik 



ni 



nj=ife!)(n-E=i^.)! 



g{si + So, S2) ^(ss + So, S2) g{si + S3, S2), 

(C.6) 

(C.7) 



denotes the multinomial coefficient. A relabeling of the dummy variables yields 

f^^\N)^y2\ ^ ^(Sl +S2,So) ^(s2 + S3,So) ^(s3 + Si,So). 

ySo, si, S2, S3 / 



From the definition (26) one easily get 



9{s,t) 



1 / n 

2 [s,t 



-1 



nA\ in A 



+ 



Therefore, ior tia — nA — n/2 we finally obtain 



n 



-So, '^l, S2, S3 



n/2 
So 



n 

l<i<3 



n 

So, Si + Sj+1 



n/2 



Now, by using the Stirling approximation and scaling the variables 



So Si S2 S3 

(70 = — , ai — — , (72 = — , (73 = — 

n n n n 



we obtain, after some algebra, the asymptotic form 

f^^\N) ~ ^2nn)^ ^ exp{n,5((7o, (Ji, (72, (73)}. 



(C.8) 



(C.9) 



(C.IO) 



(C.ll) 



(C.12) 



In Eq. (C.12) we have set (with the imphcit convention that the indices are cyclical) 

A = 



11^=1(1 -CrO-CTi- (7i+i) 



(7o(l - 2(7o)^(l - (To - (7i - (72 - (73) nj=i " 2(7^ - 2(7i+i) 



(C.13) 



and 



S'((To, (Ti, (72, (73) = 5'4((To,0'l,C72,a3) - S'2((7o,(Ti + (72) - 5'2((7o,0'2 + C^s) 



- S2{ao, (73 + (7i) + ^-Si(2(7o) + ^-Si(2(7i + 2(72) 
+ ^-51(2(72 + 2(73) + \Si{2a3 + 2(7i), 



(C.14) 



with 



n n n 

i=l 1=1 i=l 



By noting that 



dS„ 



Sn{xi, ...,Xn) = E^*^ ~ - Xj) , 



(C.16) 



1=1 



i=l 
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one easily gets 

^-oll + E-.i+^o. (C.17) 

" i=l ' 

with 



2 ^ ^i-ao-ai-a2- (73)2(1 - 2(7o)3 n?=i(l - 2(7, - 2(7i+i) 
In the hmit n — > 00 the main contribution comes from the saddle point 
(ctq, ct^, (72, CTg), solution to the system 
(99 

— = 0, with i = 0,1,2,3, (C.19) 

(9(7i 

that reads 

3 

(1 - (7o - (7i - (72 - (73) (1 - 2(7o)^ = 8(7o ]^(l - (7o - (7, - (7i+l), (C.20) 

i=l 

(1 - (7o - (7i - (72 - (73) (1 - 2(7i - 2(7i+i) (1 - 2(7i - 2(Tj+2) 

= 4(7i(l - (7o - (7, - (7i+l)(l - (7o - (7i - (7^+2), (C.21) 

with i — 1,2, 3. In the limit n — > oo we get 



,j=0 ^^'^^^ 



\d(Tido-jJ 



= A* det — ^ exp (n-S*) (C.22) 



where the starred functions A*, S*, and d'^S* /daidaj are evaluated at the saddle point 
{a*o,al,ala*s). 

The symmetry of the equations suggests to look at a symmetric solution of (C.21) 

with 

(Tj = (7 with i = 1, 2, 3, (C.23) 

which yields 

(1 - ao - 3a){l - 2aof = 8ao(l - - 2o^)^ 

(1 - (To - 3ct)(1 - Aaf = 4ct(1 - ctq - 2af. (C.24) 

We get 



^' (7^(7*3(1- 2(7^)3(1 -4(7*)3(1- (7^-3(7*) y 8(7^^*^(1 - 4(7*)^ ' ^ ^ ^ 

5* = ^ W (l-"o-2^*)^ = 1 W (1 - ^^0)' f C 26) 

2 ^ (1 - (7o* - 3(7*)2(1 - 2(7^)3(1 - 4(7*)3 2 ^ 64(7o*2 ( 1 - 4(7* ) 3 ' 

and 

/ d^S* \ _ {l-a^-5a* + 2a*{Aa* + a*)f 

\daidaj ) ~ a^a*%l - 2a*o){l - 4a*)3(l - a*^ - 3a*)(l - (t*q - 2(T*f 

X (1 - 2a* - (T*2 - 4cT*2 + 4(tV*(8(t* + 2(t* - 3)) (C.27) 
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The solution of the system that gives the largest contribution is 

^ - 197- 18v^~ 0.108955, 



36 36\/l97 - ISv^ ^6 



5043923 158161v^ 



+ 



+ 



144 (197 - ISv^) 48 (197 - 18v^) 

980473 _ 2561v^ 

72 (197 - 18v^)^ 2 (197 - 18v^)^ 

18119 563a/iT3 

+ 



144 (197 - 18\/Tl3) 48 (l97 - 18\/n3) 
~ 0.104767. (C.28) 

Plugging these results into Eq. (C.22) we get 

f^°\N) ^ cN^-^ (C.29) 

where 



c=A* det 



daidaj 

*\3 



(I -a*- 2a *) 



X 



(1 - (7$ - 5(7* + 2(7*(4(T* + a^))^Jl - 2a* - a^"^ - 4(7*2 + Aa*a*o{8a* + 2a*Q - 3) 
~ 1.05385, (C.30) 

and 



7 = 5 - ^0 log2 log2 



°^ 4.1583. (C.31) 



64af{l - 4(7*)3 



Appendix D. 

We evaluate here the asymptotic behavior of the function ji^\N) defined in (118). By 
using the definition (25) we have 

f3\N)= Yl 9{Kk2®h) g{k2,ki®h) g{h,ki®k2) 

— ^fciA(fe2efc3),0'^A:2A(fcieA:3),0<^fc3A(fciefc2),0 

X g{\ki\, \k2 e ksl) g{\k2\, \ki ® k^l) gQk^l, \ki © ks]) 
= ^ h{so, Si, S2, S3)g{si, S2 + S3 - 2so) g{s2, Si + S2 - 2so) 

X g{ss,si + S2-2so), (D.l) 

where 

h{so, 81,82,83) = ^ <^si,|fci|(^S2,|A;2|<^S3,|fe3|<^so,|feiAA;2|'^so,|feiAfe3|<^so,|A;2Afc3| 
kiMM 
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(D.2) 

and we have used 

h A {k2 ® ks) =0, k2A {ki ® /cg) =0, ks A {ki ® /cg) = 

fci A ^2 = /ci A /C3 = /C2 A ks, 

\ki ® kj\ = \ki\ + \kj\ - 2\ki A kj\ Vi,j = 1,2,3. (D.3) 

We find 

h{so, si, S2, ss) = — TT ^^77 .. . (D.4) 

Using the substitution si — >■ si + S2 — >■ S2 + Sq and S3 — >■ S3 + Sq in Eq. (D.l) we 
obtain 

f3^\^) = y2\ ^ 5'(S1 + So,S2 + S3) 5(S2 + So,Sl + S3) ^(S3 + So,Sl + S2), 

\ So, Si, S2, S3 / 

(D.5) 

in terms of the multinomial coefficient (C.7). Using (C.9 ) for ua — = n/2 we finally 
obtain 

^ ^ ySo, Si, S2, S3 y ^^.^3 + Si, Si+i + Si+2 J \So + Si J ^Sj + Si+i J 

(D.6) 

Using the Stirling approximation and Eq. (C.ll) we get 

/3^^(^) ~ ^ .2 '^('^O' '^i' ^3) exp{nS{aQ, ai, (72, (73)}, (D.7) 



where 

^crocricr2cr3ni<j<3(l - 2cro - 2ai){l - 2crj+i - 2crj+2) 



V (1 - o-Q - di - 0-2 - 0-3) 

A((7o, (7i, (72, (73) = , (D.8) 



and 

S'((7o, (7i, (72, (73) = 5'4((7o, (7i, (72, (73) — 5'2((7o + (7i, (72 + (73) 

— 'S'2((7o + (72, (73 + (7i) — 5'2((7o + (T3, "^l + ^2) 

+ ^5l(2(7o + 2(7i) + ^>5i(2(7o + 2(72) + ^^l(2(7o + 2(73) 

+ ls,{2a, + 2(72) + ^-51(2(72 + 2(73) + ^-51(2(73 + 2(7i), (D.9) 
with the entropies defined in Eq. (C.15). By (C.16) one gets (C.17) with 

log (1 - (7o - (7i - (72 - (73)^ ^ 

2 ni<i<3(l ~ 2(7o - 2(7i)(l - 2(7i+i - 2(7i+2) 

In the limit n — )■ +00 we can use the saddle point approximation. The saddle point 
is solution to the set of equations 
dS 

— = 0, with i = 0, 1,2,3, (D.ll) 
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that reads 

(1 - 2ai - 2(7i+i)(l - 2ai - 2(7,+2)(l - 2(7^ - 2ai+s) = 8(7^(1 - (Jq - (7i - ^2 - asf. 

(D.12) 

The symmetric solution, cTj = a for aU i, corresponds to the largest contribution and is 
given by 

a* = ^. (D.13) 
As in (C.22), in the hmit n — > oo we get 

f^'\N) A* det (^^) ^ exp (n^*) , (D.14) 

with 

A* = — — - = 324, det ^^^^ = — -— = 324^, (D.15) 

and 

We finally get 

/i^^(Ar) ~e"^°s(3/2) ^^a^ f^j. ^^oo^ (D_17) 

with 

a = log2 3 - 1 ~ 0.584963. (D.18) 
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